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The Classical and Quantum Mechanics of a Particle on a Knot 
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A free particle is constrained to move on a knot obtained by winding around a putative torus. 
The classical equations of motion for this system are solved in a closed form. The exact energy 
eigenspectrum, in the thin torus limit, is obtained by mapping the time-independent Schrodinger 
equation to the Mathieu equation. In the general case, the eigenvalue problem is described by 
the Hill equation. Finite-thickness corrections are incorporated perturbatively by truncating the 
Hill equation. Comparisons and contrasts between this problem and the well-studied problem of a 
particle on a circle (planar rigid rotor) are performed throughout. 
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INTRODUCTION 

The example of a particle constrained to move along a circle - the so-called planar rigid rotor - is one of the simplest 
problems that is discussed in text-books of quantum mechanics. The beguiling simplicity of this problem is at the 
heart of many non-trivial ideas that pervade modern physics. For understanding many issues like the existence of 
inequivalent quantizations of a given classical system [l| , the role of topology in the definition of the vacuum state in 
gauge theories 0 , band structure of solids Q , generalised spin and statistics of the anyonic type 0 , and the study of 
mathematically interesting algebras of quantum observables on spaces with non-trivial topology y , the problem of a 
particle on a circle serves as a toy model. 

In this paper, we consider the problem of a particle constrained to move on a torus knot. Besides adding a new 
twist to the aforementioned problems, the present system can be thought of as a double-rotor (analogous to the 
double-pendulum, but without the gravitational field) which is a genuine non-planar generalization of the planar 
rotor. 

The paper is organised as follows: In the next section we introduce toroidal coordinates in terms of which the 
constraints which restrict the motion of the particle to the torus knot are most naturally incorporated. As a warm-up, 
we then analyse the particle on a circle in toroidal coordinates. This prelude allows us to compare and contrast the 
results of the subsequent sections with the well-known results for the particle on a circle. The following two sections 
deal with the classical and quantum mechanics of a particle on a torus knot. In the penultimate section we briefly 
touch upon the possibility of inequivalent quantizations of the particle on a knot. These will be labelled by two 
parameters, in contrast to the particle on a circle. The concluding section summarises and presents an outlook. 


TOROIDAL COORDINATES 

The toroidal coordinates Q are denoted byO<77<oo, — t: < 6 < tt, Q < (j) < 27r. Given a toroidal surface 
of major radius R and minor radius d, we introduce a dimensional parameter a, defined by — d^, and a 

dimensionless parameter rjo, defined by rjo = cosh~^ (R/d). The equation rj = constant, say rjo, defines a toroidal 
surface. The combination R/d is called the aspect ratio. Clearly, larger rjo corresponds to smaller thickness of the 
torus. In the limit 770 —>■ 00 , the torus degenerates into a limit circle. 

The toroidal coordinates are related to the usual Cartesian coordinates by the equations 

asinhTycos^ asinh? 7 sin 0 asind 

^ ^ ^ y — _^ ^ ^ ( 1 ) 

(cosh ?7 — cos6) ’ (coshry — cos0) ’ (coshry — cos0) 
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The metric coefficients are given by the equations 


hi — ^2 


a 

(coshry — cos0) ’ 


/13 = /iisinh?7 


and the volume element is 


dV 


a^sinh?7 

(coshry — cosd)^ 


( 2 ) 

( 3 ) 


With the help of these basic relations, it is straightforward to rewrite well-known Cartesian expressions in toroidal 
coordinates. 


A Particle Constrained to Move on a Circle 

The Lagrangian for a free particle of mass m in Cartesian coordinates (x, y, z) is 

L=^ix^+f + z^) ( 4 ) 

In the above expression, and henceforth, an overdot refers to a time derivative. After some algebra, this expression 
can be rewritten in toroidal coordinates as 


^ m 2 + sinh^ry 0^) 

(cosh„-co.0P 

To restrict the motion of the particle to lie on a circle in the xy plane, we impose the constraints ry = ryo, a constant, 
and 9 = 9 q, another constant. The Lagrangian then takes the form 

^ ma^ sinh^ryo^^ , . 

2 (coshryo — cos0o)^ 


The Euler-Lagrange equation 


— (—) - — 
dt d<j) d(j) 


( 7 ) 


then yields, as expected. 


^ = 0 = (jjt + (j>o 


( 8 ) 


where a; is a constant and has the physical interpretation of frequency, and (j)o is a constant of integration which 
specifies the position of the particle on the circle at time t = 0 - similar to plane polar coordinates. 

Defining a rescaled mass M = m ^2 > Hamiltonian H = 2Ma^ with the momentum canonically 

conjugate to </> being given by p^j, = Ma^(j) as usual. Using this to set up the Schrodinger equation and solving it, we 
get, for the eigenvalues and the normalised eigenfunctions respectively. 


En — 


2Ma2’ 


^n(0) = 






n = 0,1,2, 


( 9 ) 


For large tjq, the thickness of the putative torus decreases and M ^ m: we approach the well-known expressions in 
plane polar coordinates. 

Interestingly, it is also possible to get a particle to move on a circle by imposing the constraints ?y = ryo, a constant, 
and (f> = (j)o, another constant. This however results in a more complicated Lagrangian viz. 


L = 


ma^ 


02 


2 (coshryo — cosd)^ 


( 10 ) 
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The resulting Euler-Lagrange equation is 

0 (cosh? 7 o — COS0) = —sin 0 0 ^ ( 11 ) 


which can be re-written as 


and readily integrated to yield 


— [ 0 (cosh? 7 o — COS0)] = 0 


0(cosh?7o — cos 9 ) = K 


( 12 ) 


K being a constant. Thus the solution is reduced to quadratures. Thanks to the presence of the factor (coshr;—cos 0 ), the 
solution is not as simple as the one in plane polar coordinates. The Hamiltonian can be obtained in a straightforward 
manner and is given by 


H = (coshryo — cosOy ( 13 ) 

The presence of the 0 -dependent multiplicative factor is portentous of additional complications that arise when we 
make a transition to quantum mechanics. In particular, the fact that the conjugate operators and 0 do not commute 
requires us to perform an operator-ordering of the classical Hamiltonian. 

The above analysis shows that while toroidal coordinates are ideally suited to consider the motion of a particle on 
a circle in the xp-plane, they are more cumbersome when it comes to handling paths which stray from the xy-plane. 
Since a knot is intrinsically non-planar, we should be prepared to confront the attendant complications. It should be 
mentioned, however, that these complications would also be present in other coordinate systems. We choose to work 
with toroidal coordinates because of their suitability in imposing the constraints that define a torus knot. 


CLASSICAL MECHANICS OF A PARTICLE ON A KNOT 

As already mentioned, the constraint ij = tjq defines a toroidal surface. A (p, q) torus knot can be obtained by 
considering a closed path that loops p times around one of the cycles of a torus while looping around the other 
cycle q times, p, q being relatively prime integers. The desired property can be enforced by imposing the constraint: 
p 9 + qcj) = 0 . It is easy to check that 0 —>■ 0 -|- 2 Trq ^ — 27 rp i.e. as we complete q cycles in the 0 direction, we 

are forced to complete p cycles in the (p direction - as required. Imposing the above two constraints on equation ( 5 ), 
we get the Lagrangian for a particle on a torus knot to be 

L=ff{cP)<P^ ( 14 ) 

where 

= l - -—M = m{a^+l 3 ‘^) ( 15 ) 

(7 — COSQ!0)^ 

with 

a = —q/p, P = sinhryo, 7 = coshryo ( 16 ) 

The main difference between the Lagrangian in ( 14 ) and the one for a particle on a circle viz. equation (6), lies in the 
appearance of the (/(-dependent factor /((/>) in the Lagrangian which contains the information about the non-trivial 
embedding of the knot in three dimensions. 

The Euler-Lagrange equation is given by 




( 17 ) 


where the prime denotes a derivative of the function / with respect to its argument (j>. Now, using the above equation 
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of motion, it is straightforward to show that 

= = 'i)=- cosacj)) (18) 

where ^ is a constant. Noting that (1 — 7 ^) < 0, the latter equation can be integrated to get 

In the limit rjQ —>■ 00 , the right hand side is linear in t, as expected for a particle on a circle. 

The momentum canonically conjugate to (j) and the Hamiltonian H can be easily worked out and are given by 
the following expressions 


p, = M/(«0, H = ^ ,20) 

QUANTUM MECHANICS OF A PARTICLE ON A KNOT 

In principle, once a Hamiltonian is given, it is a straightforward exercise to write down the Schrodinger equation. In 
the present case, the classical Hamiltonian involves a term which mixes the coordinate and the canonically conjugate 
momentum. Since these canonical pairs will be elevated to the level of operators in the quantum theory, we need 
to prescribe an ordering for the operator products. We choose the so-called Weyl ordering which symmetrises the 
product as follows: 

H =-^IjpI+p^jp^+pIj] ( 21 ) 

In the above equation, and in what follows, we refrain from putting hats, but it should be remembered that both (j) 
and p 0 are operators which obey the canonical commutation relations viz. = ih. Pulling all the momentum 

terms to the extreme right in preparation to make them act on a wavefunction, we get 

In the above form, the Hamiltonian is tailor-made for constructing the time-independent Schrodinger equation, 
with the usual prescription for replacing the momentum by the corresponding differential operator. The resulting 
Schrodinger equation is 


1 cP tP 1 ! d fP A „ 


(23) 


The first derivative in (j) can be eliminated by the well-known trick of substituting ip = in the above equation and 
getting rid of terms proportional to dY,/d4> by choosing x appropriately. This yields for x, 


X oc \/7 


(24) 


Substituting this in (23) and collecting the remaining terms, the Schrodinger equation reduces to the following equation 
for E 



+ V{(j))]T. = 0 


(25) 


where the ‘potential’ V is defined by 


V(A 


Af f-f 


'2 


-b 


2MEf 

d2 J 


12/2 


(26) 
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V is an even function of (j). Substituting for / from (15), we get after some algebra, 

2MEa‘^/h? + o? 12 — a^7COSQ;())/3 — a^cos2Q!())/6 


V = 


(7 — cosac/iy 


(27) 


Since V((/)) is a periodic function, the above potential can be expanded in a Fourier series and equation (25) gets 
identified with the Hill differential equation Q. 


The Thin-Torus Approximation 


As already mentioned, large values of tjq and hence large values of 7 , correspond to a thin torus around which the 
particle’s trajectory winds. In this limit, we can restrict to terms of the order of I/ 7 . Then /3^ "^ 7 ^, hence M —>■ m, 
and equation (27) simplifies to 


c ^2 q ,2 

V = —A — —cosad) 
4 37 ^ 


where 


Equation (25) now takes the form 


A = 


8mEa^ 


(P o? 

+ T-' - =» 

Changing variables such that a(f> = 2z, the above equation becomes 

P 4 

[ 7-7 + A - —cos22;]E = 0 
az^ 37 


(28) 


(29) 


(30) 


(31) 


which is immediately recognised to be the well-known Mathieu equation 

The solutions E of the Mathieu equation with the required periodicity are given by the Mathieu functions of 
fractional order v viz. 


, , 1 rCOs(i^ -I- 2)2 cos(^ —2)zn 

ce^{z,j) = cosvz - 


67 (ly + 1 ) 


(ly-l) 


(32) 


1 |-sin(i^-f 2)z sin(j2 — 2)z.| 

sei. z, 7 ) = sim/z - — ———- - -— 

67 '- (iy + 1) \y-^) 

The complete solution with two arbitrary coefficients A and B is given by 

E = Asei.(z, 7 ) -I- Bce^{z,'j) 


(33) 


(34) 


Setting V = ^ where n is an integer, we see that the above functions have a periodicity qn in z, which translates 
into the required periodicity 2pTT in </>. 

The condition relating A to is given by 


X = P + 


9^2(j^2_l) 


(35) 


Since we are restricting our attention to I /7 order, the boxed terms can be neglected. The allowed values of A follow 
by setting v = These values of A, in conjunction with equation (29), determine the energy eigenvalues to be 


En. = 


n ha 


2ma'^q'^ 


(36) 
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Before proceeding further, it is worth recalling that the complete solution of equation (23) that we are trying to solve 
is given by "0 = yS with y oc i//- The complete solutions for the (unnormalised) eigenfunctions 0 with the correct 
boundary conditions are therefore given by 




(7 — cosa 0 ) 


i’-'’ W = 


(7 — sina 0 ) 


X {cos{na(j)/q) 


X {sin(na0/(3') 



(37) 


(38) 


where we have used 2z = a4>. Further, since we retain only terms of order I/ 7 , the boxed terms in equations (37) and 
(38) can be neglected. 

In passing, we mention that the two independent solutions (32) and (33) can be combined into a single equation 
given by Q 


E = 




U 


(39) 


where 


u = sin( 2 : — a) + a3C0s(3z — a) + &3sin(3z — cr) + a5COs(5z — a) + 65sin(52; — u) + • • • (40) 

where cr is a new parameter such that a = ii 12 yields the solution (32) and cr = 0 yields the solution (33). In the above 
the coefficients a,b are determined in terms of 7 and cr. To order I/7 that we are interested in, only 63 = — is 
non-zero. This succinct way of writing the general solution will be particularly useful in incorporating finite-thickness 
corrections. 

It may be noted that for 9 = I, p = — 1 , and hence cr = — I, the above results (35-38) reduce to the well-known 
results for a particle on a circle. For a (2,3) torus knot, namely, the trefoil, p = 2, and the eigenfunctions have a 
period drr. The general solutions of Mathieu equations with a period 47 r were first worked out by Lars Onsager in 
1935 in his dissertation for a PhD at Yale [^. While it is gratifying to note this, it is also a little disappointing. The 
energy levels and energy eigenfunctions are the same as that of a particle on a circle, except for the factor of a. This 
is a consequence of the fact that, in the weak coupling limit (large 7), the putative torus degenerates into a limit 
circle, with the attendant vagueness associated with the winding in the 9 direction. Correspondingly, the Mathieu 
functions degenerate into trigonometric functions. It may be tempting to think that the general solution (for arbitrary 
7) will be given by Mathieu functions, with the boxed expressions in equations (35) and (37-38) being the next order 
corrections. The story, however, is slightly more complicated. Our penchant for boxing negligible pieces relates to 
this fact. 


The Slightly-Thick-Torus Correction 


To the next order in I/ 7 , the correct expression for the potential is obtained by starting with equation (27), making 
a binomial expansion of the denominator, and collecting terms up to order 1 / 7 ^. This strai ghtf orward exercise, 
followed by the steps that lead up to equation (31), yields the so-called Hill-Whittaker equation [l0| 


[-j-^ + 00 - 1 - 20icos2z + 202 COs4z]E = 0 


with 


2 

37^ 


00 — A + 




where now 


8MEa^ 

fl2^2ry2 


(41) 


(42) 


( 43 ) 
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Following Ince [1^ , the most general solution of the Hill-Whittaker equation can be obtained along the same lines 
adopted for solving Mathieu’s equation and yields the following energy eigenvalues 


En — 


%ME?- 


{-4^' + (1 - |^cos 2 u - ^sin 2 a - ^)} 


37 '’ 


97 ^ 


(44) 


The corresponding solutions are 


y;(n) _ e^™^/' 5 {sin( 2 ; — cr) -I- ^sin(32: — cr) -b — ^^)sin(5z — a) + ^cos(32; — a) — ^^cos(5z — cr)} (45) 

3 108 939 

where a is the parameter introduced earlier. Once again multiplying by the factor y; = N^/J, expanding the de¬ 
nominator, retaining terms up to order 1 / 7 ^ and, rewriting everything in terms of (j) using acj) = 2z, gives the final 
expression for the eigenstate to be 




2 . 3a0 , .iad , 2 5a^ , .bad) ,,, 

+ - cr) -b cos(— - cr) - -cosa(j){sm{— - cr) - cos(— - cr)}] 


1, . .acj) , 1 . .bacj) , 

— sin-cr -I-sm(-cr) 

7 ^ ^ 2 108 ^ 2 ' 

Acos^acj) r . .bacj) . .bacj) ,,, 
-^{sm(^ - a) + cos(^ - rr)}] 


2cosad) r . / ^ctd> X / Sad* 
^-^{sin(^ -a)+ cos(^ - a)} 


1 r .ad) . 1 . .bacj) ,, 

— [C0SQ;(()|sin(---cr) -b 77^sm(—- ct)| 


7 . 


108 


2cos^aib, . ,3ad) , .iacj) 


4cos^ad), . .bacj) , .bacj) x,,-. 
-^ - a) + cos(^ - u)}]} 


cr)} 


(46) 


where iV is a normalization constant. 


The Result For An Arbitrarily Thick Torus 

For the sake of completion, we mention that this method can be systematically continued to arbitrary orders in 
1 / 7 . The corresponding equation satisfied by E is the Hill equation given by 




4-00 + 2 ^ 02rCOs2rz]S = 0 


(47) 


r—1 


As in the earlier section, we follow Ince [l0[ , and try a general solution of the form 

E = E''^u (48) 

where 

V = pi(cr)0i +P2(cr)02 4-bqi(cr)0i +q2(cr)02 4-b ^120102 + 9130103 + 9230203 H-bri(cr)0? 4- (49) 

and 

u = sin(z - cr) + Ai(z, cr)0i + A2 (z,ct)02 • • • + Bi{z,a)Ql + B 2 {z,a)&l 4-+ Bi 2 {z, cr)0i02 4- (50) 


with cr being determined by the relation 

00 = 1 + Ai(ct)0i + A2(cr)02 • • ■ + /ri(cr)0^ + /i2(cr)02 + • ■ ■ + /il2(cr)0102 • • ■ + ICi(ct)0^ • • • 


(51) 







Substituting these expressions in equation (47) we can solve for the coefficients to any desired order, and hence obtain 
the corresponding eigenvalues and eigenvectors. We don’t pursue this exercise further since it does not shed any 
further light on the solution to the problem. 


INEQUIVALENT QUANTIZATIONS 

Let us briefly recapitulate the interesting consequences that arise if the particle which is constrained to move on a 
circle is charged, and if the circle encloses an infinitely long, infinitesimally thin, and impenetrable solenoid carrying a 
uniform current. As is well-known, the wavefunction of the particle picks up a nontrivial phase factor which depends 
on the net flux enclosed by the trajectory of the particle as it goes around the circle. Thus the wavefunction is 
multi-valued, which is a manifestation of the nontrivial topology of the circle which, in turn, is a consequence of 
the fact that the path cannot be shrunk to a point in the presence of the impenetrable solenoid. Redefining the 
wavefunction such that it is single-valued modifies the Hamiltonian in such a way that the energy spectrum depends 
on the enclosed flux. Given that the corresponding Lagrangians, with and without the flux, differ only by a total 
derivative term, the classical theory is unaltered; although different values of the flux yield different energy spectra, 
and hence inequivalent quantum theories. It is reasonable to expect similar features in the case of a charged particle 
constrained to move along a knot. 

For the torus knot of interest, two independent magnetic fluxes can be introduced. The first is the usual magnetic 
field obtained by placing a uniform current carrying, long, thin solenoid parallel to the z-axis and passing through the 
centre of the putative torus around which the knot winds. Let us denote the corresponding flux by <I> 5 . The second 
flux is obtained by a uniform poloidal current winding around the torus which produces a magnetic field which has a 
support only inside the torus, the so-called toroidal magnetic field. Let us denote this flux by $t. 

A particle constrained to move on a (p, q) torus knot, starts at a point on the surface of the putative torus and 
returns to the initial point after completing one circuit of the knot; in the process winding around the solenoidal flux 
p times and the toroidal flux q times. The total flux enclosed is therefore: $ = -|- q^T- Thus we have the equation 

which highlights the multi-valued nature of the wavefunction viz. 

0 + 2Trq, (j) - 2ttp) = exp(i$)V’(? 7 o, 0, (j)) (52) 

Defining the single-valued wavefunction 

$ 

V’(??o, 6 *, 4>) = exp(-i-—<(’) (53) 

Ipir 

and the corresponding Hamiltonian obtained by the transformation 

H = (54) 

2p7T 2p7T 

we see that the momentum operator in the Hamiltonian is shifted by which leads to a corresponding shift in 
v and hence the energy spectrum defined in equations (36) and (44). It is noteworthy that the phase picked up by 
the wavefunction of the particle, for a given (p, q) knot, is a sum of two independent fluxes. Thus the inequivalent 
quantizations are labelled by two parameters. 


CONCLUSIONS AND OUTLOOK 

The classical and quantum mechanics of a particle constrained to move on a torus knot were studied. The results 
were compared and contrasted with the well-known results for a particle constrained to move on a circle. Defining 
the knot as a trajectory which winds around a putative torus in a well-defined fashion, and using toroidal coordinates 
to parametrise the knot, makes it possible to rewrite the time-independent Schrodinger equation as a Hill equation 
which can then be studied perturbatively in the thickness of the putative torus. 

Attributing a charge to the particle and introducing two independent magnetic fields having supports in physically 
disconnected, but topologically linked, regions, leads to a two-parameter family of inequivalent quantizations of the 
particle moving on a knot. 
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The model discussed in this paper has several features which are worth discussing further. First, it would be 
natural to study the model non-perturbatively i.e. using instanton methods made popular in ii- Second, it would 
be interesting to generalise the treatment to more than one particle moving on the knot. The non-trivial phase 
factor can then be related to exotic quantum statistics of the anyonic type. It would also be interesting to construct 
coherent states and study algebras of quantum observables associated with a particle on a knot. All these problems 
have natural analogues for the corresponding, but much simpler, example of a particle constrained to move on a circle. 

I thank G. Krishnaswami and A. Laddha for discussions. This work is partially funded by a grant from Infosys 
Foundation. 
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